Stability of shortest paths in complex networks with random edge weights 
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We study shortest paths and spanning trees of complex networks with random edge weights. Edges 
which do not belong to the spanning tree are inactive in a transport process within the network. The 
introduction of quenched disorder modifies the spanning tree such that some edges are activated and 
the network diameter is increased. With analytic random- walk mappings and numerical analysis, 
we find that the spanning tree is unstable to the introduction of disorder and displays a phase- 
transition- like behavior at zero disorder strength e — 0. In the infinite network-size limit (N —* 00), 
we obtain a continuous transition with the density of activated edges $ growing like $ ~ e 1 and 
with the diameter-expansion coefficient T growing like T ~ e 2 in the regular network, and first- 
order transitions with discontinuous jumps in $ and T at e = for the small-world (SW) network 
and the Barabasi-Albert scale-free (SF) network. The asymptotic scaling behavior sets in when 
N S> N c , where the crossover size scales as JV C ~ e~ 2 for the regular network, N c ~ exp[ae~ 2 ] for 
the SW network, and N c ~ exp[a| lne|£~ 2 ] for the SF network. In a transient regime with N -C N c , 
there is an infinite-order transition with <E> ~ T ~ exp[— a/(e 2 In N)] for the SW network and 
~ exp[— a/(e 2 In N/ In In N)] for the SF network. It shows that the transport pattern is practically 
most stable in the SF network. 

PACS numbers: 89.75.-k, 05.10.-a, 75.10.Nr, 05.40.Fb 



I. INTRODUCTION 



A network is a new paradigm to study complex sys- 
tems in many disciplines in science [0,0]. A complex 
system consists of a large number of interacting units, 
and the nature of the interaction determines equilibrium 
and dynamical properties of the system. Frequently the 
simplifying assumption is made that the units are ar- 
ranged to form a simple pattern like a regular lattice or 
to interact with all others as in a mean-field theory. Re- 
cent studies, however, have revealed that the structure 
of complex systems is much richer |2j |3j. In general 
this structure is captured by a network which consists 
of vertices representing the units and edges connecting 
interacting vertex pairs. 

Complex networks exhibit so-called small-world phe- 
nomena: vertices are highly clustered and the average 
separation between vertices grows slowly with the total 
number of vertices. Watts and Strogatz Q introduced 
a small-world (SW) network as a model for these phe- 
nomena. It is obtained from a regular lattice with edges 
randomly rewired with probability p r . Later it was found 
that some complex networks have a power-law distribu- 
tion Pdeg.{z) ~ z~ 7 of the degree z. The degree of a 
vertex is the number edges incident upon it. The class of 
networks with a power-law degree distribution is called 
the class of scale- free (SF) networks and is found in many 
areas including physics, computer science, biology, soci- 
ology, etc (we refer readers to Ref. for examples) . The 
Barabasi-Albert model Q generate a SF network (7 = 3) 
growing via a preferential attachment rule H. Initially 
one starts with Zq vertices and introduces a new vertex 
at each step [growing]. It is then attached to Z exist- 
ing vertices, which are selected with probability linearly 
proportional to their degree [preferential attachment] . 



The discovery of the new classes of networks trig- 
gers extensive research. Order-disorder phase transi- 
tions H [?J U and nonequilibrium phase transitions |J 
have been studied. Interestingly, the critical behavior 
is described well by mean-field theories and strongly de- 
pends on the degree distribution R] . Stability of complex 
networks has also been studied against a strong disorder 
such as a vertex dilution [ll], [l^] . When the fraction 
of diluted vertices increases, a network may disintegrate 
into finite clusters undergoing a percolation-type transi- 
tion. In this paper, we study the effect of weak disorder 
on the transport properties of networks. 

The shortest path plays an important role for the trans- 
port within a network Jl3| , p"4| . A path denotes a sequence 
of vertices, successive pairs of which are connected via 
edges. In general there exist many paths connecting two 
given vertices. The shortest path is the one with min- 
imum path length among all the paths. The minimum 
path length is called a separation between the two ver- 
tices. Suppose one needs to send, e.g., a data packet 
from one computer to the other through the Internet. 
The shortest path provides an optimal path way since 
one would achieve a fastest transfer and save system re- 
sources. The shortest path is also important in studying 
an internal structure of a network [lq ]. The separation 
can be used as a measure of intimacy between vertices. 
The number of shortest paths that pass throu gh a vertex 
is called the "betweenness" or "load" @ @, ||, []]J. It 
reflects the importance of a vertex in mutual relationship 
or in transport. The load follows universal power-law dis- 
tributions in scale-free networks Jl^, [lj| . 

Consider shortest paths from a vertex s, called the 
source, to all other vertices in a network. In an un- 
weighted network the length of a path is just the number 
of edges that it contains and shortest paths can simply 
be found using the "breadth-first search algorithm" |L7| . 
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FIG. 1: A network (left) and its spanning tree (right) T s . 
Vertices in T s are arranged in the hierarchical order of the 
separation /. The vertex Ve, and hence V5, has twofold de- 
generate shortest paths, so T s contains a loop. Dashed lines 
represent edges which do not belong to T s . 



A sub-network consisting of all the shortest paths from 
s is the spanning tree T s , which characterizes the opti- 
mal transport pattern. Figure |l] shows an example of a 
network and its spanning tree. It is convenient to rep- 
resent the spanning tree by a diagram in which vertices 
are arranged hierarchically in the ascending order of their 
separation from the source. Then, the shortest path to 
each vertex is given by a directed path on T s . In general, 
the spanning tree does not have a tree structure. If there 
are degenerate shortest paths, T s contains a loop. Some 
edges do not belong to T s . They do not contribute to 
any flow from/to s. 

It has been assumed that all edges are equivalent hav- 
ing the same cost. However a real network would be 
described better with weighted edges. A weight of an 
edge may represent an access cost, a p hysical length, 
or an intimacy between vertices fl5[ |l8[ . For example, 
edges between scientists in scientific collaboration net- 
works may have weights which depends on the number of 
coauthored papers p5| . In a weighted network, the path 
with the minimum number of edges is not necessarily an 
optimal one. In this work, we study disordered networks 
with randomly-weighted edges and investigate the effect 
of the disorder on the transport pattern. The regular 
network, the SW network Q, and the Barabasi- Albert 
SF network are considered. This paper is organized 
as follows: In Sec. 0, we define the shortest path and the 
spanning tree in the disordered network. The disorder 
modifies the shape of the spanning tree. The response is 
described for the regular and SW networks in Sec. Ill and 
for the SF network in Sec. |y|. We conclude in Sec. V|. 



II. SHORTEST PATHS AND THE SPANNING 
TREE OF DISORDERED NETWORKS 

Consider a disordered undirected network. An edge 
e between two vertices u and v will be denoted as e = 
[u; v] = [v; u]. To each edge e a non-negative weight c(e) 
is assigned which is called the edge cost of e. Here we 
neglect all system-dependent details and assume that 



c(e) = 1 + 57(e), 



(1) 



where 77(e) 's are random variables distributed indepen- 
dently with distribution !F(ri) (77(e) > — 1). 

In a disordered network, a minimum- cost path plays 
the role of the shortest path in a pure (77 = 0) network. 
For given vertices u and v, the minimum-cost path is 
given by the one with minimum path cost. The path cost 
is defined as the sum of all edge costs in the path. With- 
out disorder (77(e) = for all e), the path cost is equiva- 
lent to the path length and the minimum-cost path is the 
same as the shortest path. Hereafter, the minimum-cost 
path will be called the shortest path, and the minimum 
cost is denoted as a distance. The path length of the 
shortest path is called the separation of the two vertices 
connected by it. 

The spanning tree T s of a disordered network can be 
found using the "Dijkstra algorithm" jL7| : Divide all ver- 
tices into two sets S and its complement S. Initially 
S = {s} and the source is assigned to a distance label 
d(s) = and a separation label l(s) — 0. At each it- 
eration, one selects an optimal edge e* = [u*;v*] that 
has a minimum value of d(u) + c([u; v]) among all edges 
e = [u; v] with u € S and v € S. Then, the vertex v* gets 
the labels d{v*) = d(u*) + c{[u* ; v*]) and l(v*) =l(u*) + l 
and a predecessor label pred(v*) — u* , and is shifted 
from S to S. The iteration terminates when the set S is 
empty. The shortest path to each vertex is then found 
by tracing the predecessor iteratively back to s. The dis- 
tance and separation from s to each vertex v are given 
by d(y) and l(v), respectively. The average separation 
^ >s = Iv^T Eu^s l( v ) wm " be called a diameter. 

In a homogeneous (e.g. regular) network and a weakly 
disordered network (e.g. SW network), all vertices are 
equivalent after an average over all disorder realizations. 
The diameter is independent of the source, D s = D. In 
such cases we select the source s arbitrarily. On the other 
hand, the SF network has a highly inhomogeneous struc- 
ture. For the SF network, we select the hub which has 
the largest degree as the source s since it plays the most 
important role in the transport |Dj ] . 

The spanning tree T s of a disordered network will be 
different from T° of the same network without disorder. 
If the disorder has a continuous distribution, with prob- 
ability one all shortest paths are uniquely determined. 
Therefore, T s has a tree structure, whereas T° may have 
loops. Without loops, T s consists of only (N — 1) edges. 
Moreover, a vertex may have a shortest path that can- 
not be found in T°. For example, a path (s, V\, vg, Ug, v@) 
in a network in Fig. |l| may be a shortest path to vq, so 
that T s includes the edge [vg; vq] which does not belong 
to T°. Such an edge of T s that does not belong to T° 
will be called an activated edge. The disorder activates 
it to play a role in the transport. The activated edge 
results in a drastic change in the shape of the spanning 
tree and increases the network diameter. We quantify 
the change by the density of disorder-induced activated 
edges $, which is given by the number of activated edges 
in T s divided by (TV — 1), and the diameter-expansion co- 
efficient T = (D s — D°)/D° with D s (D° s ) the diameter 



3 




(a) 



(b) 



(c) 



( d ) 



FIG. 2: A regular network with Z = 2 in (a) and its spanning 
trees without disorder in (b) and with disorder in (c) and (d). 
(c) shows a spanning tree of type B and (d) shows a spanning 
tree with both segments of type A and type B. Vertices in 
(b), (c), and (d) are arranged in the hierarchical order in the 
separation from the source s. Activated edges are represent 
by thick lines in (d) . 



with (without) disorder. The activated edge emerges as a 
result of competition between all paths connecting a ver- 
tex to the source. So networks with different structures 
respond differently. In the next section we will study the 
evolution of the spanning trees of the regular network 
and the small- world network [0 . 



III. SMALL- WORLD NETWORK 

Consider a regular network consisting of N vertices on 
a one-dimensional ring, each of which is connected up 
to Zth nearest neighbors with undirected edges. The 
small-world (SW) network is obtained by rewiring each 
edge with probability p r (see Ref. || for a detailed pro- 
cedure). Except for extreme cases with p r — (regular 
network) and p r — 1 (random network), the SW network 
displays the small- world phenomena Q. We introduce 
a quenched disorder in edge costs as in Eq.|l]) with the 
disorder distribution 







for — e < rj < e 
otherwise. 



(2) 



Disorder strength is controlled by the parameter e (< 1). 
First we focus on the regular network (p r = 0) with Z = 
2, which gives us a lot of insights. 



A. Regular network (p r — 0) with Z = 2 

A regular network with Z = 2 and N = 18 is shown in 
Fig. ||(a). We consider the shortest path from the source 
s. Without disorder (e = 0), the shortest path to b n is 
unique. On the other hand, there are n degenerate short- 
est paths to a n . So the spanning tree T° has a ladder 



shape with diagonal rungs from b n to a„+i as shown in 
Fig. |(b) g§. All edges {[a n ;b n ]} are missing in T s . 

With infinitesimal disorder, all loops break up in the 
spanning tree. So either [a„_i;a ra ] or [6 n _i;a„] should 
be removed from T s . Consequently, T s has a tree shape 
with a single branch for b n 's and with sub-branches for 
a„'s, cf., Fig. ||(c). This shape will be denoted as a type-B 
tree. The branching points are determined from recur- 
sion relations for the distances d(a n ) and d(b n ): 

d(a n ) = min{d(a„_i) + c([a„_i;a„]), 

d(b n -i) + c([6 n _i;a n ])} (3) 
d(b n ) = d{b n -i) + c{[b n -i;b n \) . (4) 

If d(a„_i) + c([a„_i;a„]) < d(b n -i) + c([6 n _i; a n }), the 
predecessor of a n is a„_i, otherwise the predecessor is 
b n -\. The recursion relation (^^) holds if T s has the 
type-B structure, which is valid as long as 



d(a n ) + c([a n ; b n }) > d{b n ). 



(5) 



[a n ; b n ] 



When it is violated, the predecessor of b n is a n , 
becomes activated, and T s changes its shape. 

We introduce now a random walk interpretation of 
the recursion relation. Define Xa(ji) = d(a n ) — n and 
X B = d(b n ) — n and insert it into Eq. (||). Using Eq. (Q), 
one obtains A"s(n) = A"s(n — 1) + r;([6„-i; 6„]); Ib(ii) 
can be interpreted as the coordinate of a one-dimensional 
random walker (walker B) after n jumps. XA(n) is 
given by the minimum of X^(n — 1) + ?7([a n _i; a n ]) and 
Xb{ti — 1) + r]([b n -i- a n ]) ■ The first term also suggests 
that XA(n) can be interpreted as the coordinate of an- 
other random walker (walker ^4) after n jumps. But its 
motion is constrained: After each jump, one has to com- 
pare the current position of A [the first term] with the 
position of B [the second term~ Xs(n) + 0(e)], and take 
the minimum as Xa{ti). Therefore, one may assume a 
hard-core repulsion between two walkers (the interaction 
range fluctuates by an amount of 0(e)). The inequal- 
ity in Eq. (|^) can be rewritten as ATs(n) — XA(n) < 
c([a„; b n ]) — l + 0(e), which imposes a constraint on the 
random-walk motion. 

It is convenient to introduce Xji{n) = Xb{tl) — 
XA{n). It can be interpreted as the coordinate of a one- 
dimensional random walker R in the presence of a fluc- 
tuating reflecting wall at Xr — 0(e) and a fluctuating 
absorbing wall at Xr = 1 + O(e). At each time step, the 
walker R performs a jump of size 77' obeying the distri- 
bution T R {rj') ee J drjif drj 2 T(i^i)T{ri2)5{rj' - rjx + r) 2 )- 
Hereafter the boundary walls are assumed to be fixed at 
Xr = and Xr — 1. The fluctuations do not modify 
the scaling behavior of $ and T with a possible change 
in the coefficient of the leading order term. 

The random walker R determines the shape of the 
spanning tree. If the walker does not touch either wall 
at Xr — or 1 at a moment n, then the predecessor of 
a„ (b n ) is a„_i (6 n _i). When it bounces at the reflect- 
ing wall in step n, a n has as its predecessor, and 
the spanning tree has a new sub-branch, cf. 03, 04, and 
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a7 in Fig. ||(c). If it collides with the absorbing wall at 
step no, the inequality (0) is violated and b no has a na 
as predecessor instead of b no -i (see the spanning tree in 
Fig. ||(d) which has n = 5). 

The same random-walk mapping can be established 
after an edge [a no ',b no ] is activated. If one interprets 



Xn i bn ] 

as a new source and re-defines a' 



and 



'n+no— 1) 



the distances from s 1 to a' and b' n satisfy 
the same recursion relation as in Eqs. (0) and (^) and 
the same constraint as in Eq. (||) with a (b) replaced by 
b' {a'). Thus, the spanning tree consists of a single branch 
for a"s and sub-branches for 6"s. This shape will be 
denoted as a type A, see Fig. ||(d). The creation of sub- 
branches and the switch into type-A tree are described 
by the same random-walk mapping used for the type-B 
segment. 

Combining the mappings for type-A and type-B seg- 
ments, the shape of the whole spanning tree can be deter- 
mined by the the random walker R in the presence of two 
hard-core walls at Xr — and 1. Initially ATr(O) = 0, 
and the wall at X R ~ (1) is reflecting (absorbing). 
When the absorbing wall is at Xr = 1 (0), the span- 
ning tree has the type-B (type-A) shape. The sub-branch 
emerges when the random walker collides with the re- 
flecting wall. When it collides with the absorbing wall, 
the role of the two walls are exchanged and the span- 
ning tree switches its shape. It is interesting to note that 
the random walk with two types of boundary walls were 
used to find the exact ground states of one-dimensional 
random-field Ising-spin chain (2(]] . 

The shape of T s is characterized by two length scales £i 
and £2, see Fig. ||. The former characterizes the length of 
the sub-branch, and the latter the length of each type-A 
or type-B segment. They are given by the mean time 
scales between successive collisions with the reflecting 
wall and the absorbing wall, respectively. Then, £ 1 can be 
approximated as the life time of the random walker R, be- 
ing at the origin initially, in the presence of two absorbing 
walls at X = and X = 1. And £2 is given by the mean 
life time of the random walker R, being at the origin ini- 
tially, in the presence of two absorbing walls at X = ±1. 
Such time scales are calculated in time-continuum limit 
in Appendix 



Using the results in Eqs. 



) and (|AJ) 
, we obtain 



with a = a and a 2 = j dvj ' Tnij]')^' 2 

that £1 ~ ^3/2 e _1 and £2 = (3/2)e~ 2 . Note that the 
0(e) fluctuation of the walls does not change the scaling 
exponents, but may modify the coefficient of £1. 

One edge among Z edges in a row is activated when 
a sub-branch appears only in the type-A segment (see 
Fig. 0(d)). So, the activated-edge density is inversely 
proportional to £1 : 
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where Z = 2 and (Z — 1)/Z is the fraction of the type-A 
segments. When the spanning tree changes its shape at 
a certain vertex (cf. a 5 in Fig. §(d)), its all descendent 
vertices increases their separation from s by one. So, the 
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FIG. 3: reg in (a) and Treg in (b) for the regular net- 
work with Z = 2 and N = 10000 (o), 20000 (rj), 40000 (©), 
and 80000 (A). The straight lines represent 0.24a; in (a) and 
0.32a; 2 in (b). Insets show Z dependence of the coefficients 
X* (&reg ^ X*e) and xt (Treg ^ Xr£ 2 )- 



diameter-expansion coefficient is inversely proportional 
to £2: 



1 



REG\£) 



1 



(7) 



The results are valid in the asymptotic limit where 
N 3> max[£i,£ 2 ] = £2, which suggests the finite- 
size-scaling form $ RE q(e,N) — N~ 1 / 2 Q$(eN 1 / 2 ) and 
T REG (e,N) = N~ 1 g- T (eN 1 / 2 ). The scaling functions 
behave as Q$(x) ~ x$ x an d Qt(x) ~ Xtx 2 for x 1. 
The scaling behavior is confirmed numerically. We com- 
pute both quantities for the regular network with Z = 2 
and N = 10000, . . . , 80000, which were averaged over 200 
samples. They are plotted in Fig. ||, where data collapse 
very well. From a least-square fitting we obtained that 
X$ — 0.24 and yt — 0.32, which is close to the analytic 
results (||) and @j. 

The extension to the regular networks with Z > 2 
is straightforward. Vertices are labelled as (s, v\ , . . . , vf , 
v\, . . . , V2 , ■ ■ •) starting from the source s. Then, without 
disorder, the spanning tree T° s has a Z-leg ladder struc- 
ture with diagonal rungs. Each leg (i = 1, . . . , Z) consists 
of vertices {v l n }. The predecessors of a node v l n+1 are v 3 n 
with j = i, . . . ,Z. Edges {[v l n \ v J n ]} with i ^ j do not 
belong to . When the disorder turns on, there emerge 
activated edges. Until one finds an activated edge, the 
distance to each vertex from the source satisfies recursion 
relations 

d(vi) = . min {d(<_ x ) + c(K_ x ;<])} . (8) 

j=i,. ..,Z 

The recursion relations are valid as long as 

d(vi)+c([vi;vi])>d(vi) 



(9) 



for all i 7^ j. With the mapping Xi(n) = d(v l n ) — n, one 
can interpret Xi(n) as a coordinate of a random walker 
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Ai after n jumps (each jump has the distribution J-(r))). 
Then, Eq. (|§[) implies a hard-core interaction between 



walkers, so Ai cannot overtake Aj's with j > 



Ef- 



fectively, it suffices to consider the hard-core interaction 
only between Ai and Ai + \. If Ai and Ai+i collide at 



step n, then the vertex v % n takes 



as its predecessor. 



Otherwise, v l n _ 1 is the predecessor of v\. The constraint 
(0) implies that the relative distance between all walkers 
should be less than f + 0(e), that is, \Xi - X z \ < 1. 
An activated edge [<x^ D ;o^ J appears when this inequal- 
ity is violated in the (no + l)th step. Then, we can use 
the same random-walk mapping after a cyclic permuta- 
tion (Ai, . . . , Az) — > (A 2 , . . . , Az, Ai), which continues 
repeatedly. As in the case with Z = 2, the shape of 
the spanning tree is characterized by the length scale 
£i, the mean time scale for a collision between adjacent 
walkers, and £2, the mean time scale for violating the 
constraint \Xi — X z \ < f . The time scales are approx- 
imately equal to those for the two-random-walker prob- 
lem with a constraint \X\ — X2 | ^ 2/Z. With this ap- 
proximation, one gets £1 ~ (Ze)^ 1 and £2 ~ (Ze)~ 2 . 
Using $reg = fs^- and T REG = ^, we finally ob- 
tain $jj£G — X* £l ancl ^REG — Xx£ 2 with Z-dependent 
coefficients x* oc (1 — ^) and xt oc Z. The scaling 
exponents are universal for all Z. We determine the co- 
efficients of e and e 2 for $reg an d Treg, respectively, 
numerically and plot them as a function of Z in the insets 
of Fig. |[ One sees that x* oc (1 — 1/2) and xt oc 2, as 
estimated above. 

We conclude that quenched disorder is a relevant per- 
turbation to the spanning tree of the regular network. 
Using the random-walk mapping, we have shown that a 
finite fraction of edges in the spanning tree are modified 
at nonzero disorder strength. This fraction is linearly 
proportional to the disorder strength, &reg ~ £■ We 
have also shown that the diameter-expansion coefficient 
is proportional to the square of the disorder strength, 
^ REG ~ e ■ 



B. Small- world network 

Next we study the effect of the quenched disorder in 
the SW network. Figure || shows an example of a SW net- 
work and its spanning tree T° without disorder (e = 0). 
The random rewiring of edges randomizes Tj, too, which 
is not suited for an exact description. Therefore we study 
&sw and Tsw first numerically. We calculate $>sw m 
SW networks with Z = 4 and p r — 0.2 and compare 
it with &reg m Fig. |[ &sw shows two noticeable fea- 
tures: (i) At small e, &sw appears to display a threshold 
behavior at a nonzero value of e. (ii) $>sw has a strong 
size dependence. $sw does not approach an asymptotic 
saturation at any value of e and N considered. The same 
feature is commonly observed at other values of Z and 
p r , and also for Ysw- I n what follows we describe the 
response of the spanning tree with an effective random 
walk process and we will show that the origin of the ap- 




FIG. 4: A small-world network (left) with N = 40, Z = 2, and 
p r =0.1 and its spanning tree T° (right) without disorder. 
The dashed lines indicate edges which do not belong to T°. 
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FIG. 5: Activated edge densities in the regular network &reg 
with Z = 4, the SW network &sw with Z = 4 and p r = 0.2, 
and the SF network 3>sf with Z = Zq = 4. 



parent threshold behavior is the presence of a well-defined 
e-dependent crossover size in the network. 

All edges that connect vertices at the same hierarchy 
level in T°, as represented by dashed lines in the ex- 
ample in Fig. are candidates for an activated edge. 
Focus on a pair of vertices a and b with [a; b] ^ T° 
in Fig. [|. They are descendants of a common ances- 
tor s'. The edge [a; b] belongs to T s if a difference be- 
tween costs of the two paths (denoted by thick lines) 
from s' to a and to b is larger than the edge cost 
c([a;b]) = 1 + rj([a;b]) = 1 + 0(e), where 0{e) term 
can be neglected. The probability with which this hap- 
pens will be denoted as P ac t.{[a; b}). Let l(a,b) be the 
separation of a and b to their common ancestor s' in 
T°. Using Eq. (Q), the path costs are given by a sum 
of Z(a, b) independent random variables rj's (plus l(a, b)). 
So, with the common term l(a,b) discarded, they can be 
interpreted as coordinates of two one-dimensional ran- 
dom walkers after l(a, b) jumps. Each jump follows the 
distribution ^{rj). P ac tX[ a 'i ^ s then given by the prob- 
ability that the distance between two random walkers 
is larger than 1 after I jumps. Or equivalently, it is 
given by the probability that a random walker devi- 
ates from a starting position by a distance larger than 
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1 after I jumps, where each jump follows a distribution 
•?W) = fdr)ij mHviMV2)S(v' -V1 + V2)- Theprob- 
ability distribution of the random walker after I steps 
is given by P R . w .(x,l) = {2i:a 2 l)- x l 2 e-^l {2 ^^ with 
a 2 = J dr]'ri' 2 !Fji(rf) — 2e 2 /3. Therefore, one obtains 



Pact, ([a; b]) = erfc 



(10) 



is the complementary error 



where erfc(a;) 
function. 

Now we make a mean-field-type approximation that 
each edge e ^ T° may be an activated edge indepen- 
dently. That is only true for the edges {e^ = [a^; bi] ^ T"} 
if the paths from ai and bi to their common ancestor 
do not overlap for different e^'s. Otherwise, two proba- 
bilities Pact.(e) and P ac t.( e ') are n °t independent. If e 
is activated, it modifies the shape of the spanning tree 
and hence P a ct.( e ')i and vice versa. The approxima- 
tion would give the correct scaling behavior if the span- 
ning tree would be random and self-averaging. In the 
mean-field scheme, the activated edge density is propor- 
tional to the probability in Eq. ([l(]) with l(a, b) replaced 
by the average value, i.e., the diameter of the network 
D„~lnN [El: 



&sw(N, e) ~ a\ erfc 



a 2 



eVhTiV 



(11) 



with constants a\ t 2 being independent of e and N . 

The scaling function behaves as erfc (a;) ~ 1 — -^=x for 
small x. So, in the asymptotic limit where N ^> N c ~ 
e" !£ , the activated edge density is a constant {§sw — 
Qi] in the infinite (N — > oo) network for all values of 
e 7^ 0. Therefore, the spanning tree of the SW network 
undergoes a discontinuous transition at e = e c — in the 
asymptotic limit, which is contrasted to the continuous 
transition in the regular network. 

The asymptotic behavior sets in when the network size 

2 — 2 

is bigger than the crossover size iV c ~ e~ a ^ e . It grows 
very fast as e goes to zero. For instance, when e = 0.1 
with a 2 ~ 1, only SW networks with N > N c ~ 10 43 are 
in the asymptotic region, which is improbable in any real 
networks. So, the behavior in the non-asymptotic regime 



is more important in practice. Using erfc(x) 
for x 3> 1, we obtain that 



,. aiS\/lnN 
®sw 1=- exp 



e 2 In AT 



(12) 



for -c N c . It has an essential singularity at e = and 
increases continuously and extremely slowly. Therefore, 
the spanning tree of finite SW networks undergoes an 
infinite-order transition at e = 0. 

Numerical data are in good agreement with the mean- 
field results. In Fig . |6|, w e plot $sw as a function of the 
scaling variable eVln N. Data at different values of 
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nJV] 1 ' 2 according to 
a + b/hxN for fixed 
values of e, following Eq. (|l^). The straight lines are obtained 
from a least-square fitting. 



FIG. 6: Scaling plot of $sw vs. e 
Eq. (|ll]). The inset showsjhat ln$ 



collapse onto a single curve, which supports the scaling 
form in Eq. ([TT]) . In the inset, we plot $>sw at given 
values of e against 1/bxN for N = 100, . . . ,800000 in 
the semilog scale. We observe that the data align along 
straight lines, which support the result in Eq. (|l2|). Ex- 
trapolating the straight lines to the limit In N —> oo, we 
obtained that $siy's converge to a\ ~ 0.25. 

In the regular network, an activated edge may or may 
not increase the separation of vertices from the source, 
which leads to the different scaling behaviors of Qreg ~ 
e and Treg ~ £ 2 - However, in the SW network, all 
activated edges considered above increase the separation. 
So T sw shows the same scaling behavior as &swi which 
we also confirmed numerically. 



IV. B ARAB ASI- ALBERT SCALE-FREE 
NETWORK 

In this section we study the response of the spanning 
tree of a SF network to the quenched disorder. We con- 
sider the Barabasi- Albert (BA) network with the degree- 
distribution exponent 7 = 3 (|. Among all vertices, the 
hub which has the l arg est degree plays an important role 
in the SF network [ pL 3| . So we consider the spanning tree 
of the BA network with the hub as the source. is 
measured for the BA network with Zq = Z = 4, and 
compared with Qreg and $sw in Fig. ^. $sf behaves 
similarly to &sw- 

Though the SF network is different from the SW net- 
work in many aspects, its spanning tree also has a ran- 
dom structure. Thus, we expect that the scaling behavior 
of &sf can be explained with the same mean-field-type 
approximation as <&sw- Under the approximation, each 
edge that does not belong to T° may be activated in- 
dependently. Then, the mean density of the activated 
edges is obtained by replacing l(a, b) in Eq. ( JlCj ) with the 
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diameter of the network. Recently, it was found that the 
scale- free network is ultra small in diameter J2S| . In par- 
ticular, for 7 — 3 (7: the degree-distribution exponent) 
the diameter scales as D ~ In N/ In In A, i.e., with a dou- 
ble logarithmic correction to the SW scaling D ~ In A. 
It leads us to conclude that 



§>SF = &i erfc 



a.-, 



^lnA/lnlnA^ 



(13) 



with constants and a^. 

This scaling form is indeed supported by our numer- 
ical data. In Fig. [?], we plot $>sf for the BA network 
with Z = Zq = 2 as a function of the scaling vari- 
able £y/\a.N/ In In AT, and obtain a good data collapse. 
Though the double-logarithmic correction is very weak, 
the data do not collapse at all without it. We note that 
the data also scale well with a scaling variable e(ln A) 1 / 4 , 
however, we presume that this is accidental since the ra- 
tio between (In N/ In In A) and (In A) 1 / 2 remains almost 
constant up to A = 800000. The inset in Fig. @ shows 
that ln$5i? is linear in l/(e 2 In A/ In In A). It implies 
that the scaling function has the essential singularity as 
111 Eq. @. The same scaling behavior is observed uni- 
versally for other values of Z = Zq = 3 and 4, and also 
for T sf- 

The SF network shares the same scaling behav- 
ior with the SW network with the scaling variable 
e[lnA/ In In TV] 1 / 2 . The spanning tree of the SF net- 
work undergoes a discontinuous transition in the asymp- 
totic regime: <&sf an d Ysf in the infinite SF networks 
jump from zero at e — to finite values as the disorder 
turns on. The asymptotic behavior sets in only when 
A ^> N c ~ exp[a| lne|e~ 2 ] with a positive constant a. 
In the transient regime with A <C N c , the spanning tree 
of the SF networks undergoes a continuous infinite-order 
transition with the essential singularity in $sf and T sf 
at s = 0. Note that the crossover size N c of the SF net- 
work grows much faster than that of the SW network. 
Therefore we conclude that the scale-free network has 
the most stable transport pattern. 



V. CONCLUSION 

In a network the shortest paths between vertices 
play an important role. The optimal transport pattern 
from/to a vertex s is characterized by the spanning tree 
T s , a set of all shortest paths. We have investigated 
the response of the spanning tree of complex networks 
to a quenched randomness in edge costs, and found that 
quenched disorder is a relevant perturbation. As the dis- 
order turns on, the spanning tree evolves with a finite 
fraction of edges modified and the network diameter ex- 
pands. 

For the regular network, the shape of the spanning tree 
can be described exactly using the random-walk map- 
ping. The spanning tree undergoes a continuous transi- 
tion with Qreg ~ £ and T ' reg ~ £ 2 with e the disorder 
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FIG. 7: Scaling plot of $ SF vs. e[ln N/ In In iV] 1/2 for 
the Barabasi- Albert networks with Z — Zo — 2 and 
N = 20000, 800000. Inset shows a plot of $ S f vs. 
l/(e 2 [ln N/ In In N]) in semilog scale. 



strength. For the SW and SF networks, we obtain the 
scaling form <&(e,A) = ^(eD 1 ' 2 ) where D is the net- 
work diameter. It scales as D ~ In A in the SW net- 
work and D ~ In A/ In In A in the BA network [§| . The 
diameter-expansion coefficient satisfies the same scaling 
form. The scaling function approaches a constant as its 
argument x = eD 1 / 2 — > 00. It has the essential singu- 
larity, <&(x) ~ exp[— as" 2 ] with a positive constant a, as 
x — > 0. Therefore, <!> and T have a discontinuous jump 
at e = in infinite-size SW and SF networks. It shows 
that the SW network and the SF network are more af- 
fected by the quenched disorder than the regular network. 
However, the asymptotic scaling behavior emerges only 
when the network size is larger than the crossover size 
A c . It grows as A c ~ exp[ae~ 2 ] for the SW network and 
A c ~ exp[a| lne|e -2 ] for the BA network with a positive 
constant a as e vanishes. Since the crossover size grows 
extremely fast as e goes to zero, the behavior in the tran- 
sient regime with A <C A c is more important in practice. 
In contrast to the discontinuous jump when N ^S> N c , $ 
at A <gC A c grows continuously and very slowly with the 
essential singularity at e — 0. Numerically we obtained 
that $ < 10" 7 for e < 0.1 up to A = 800000. Therefore, 
we conclude that the transport pattern in finite SW and 
SF networks is extremely robust against the quenched 
disorder. In particular, the SF network is most stable. 

The interesting scaling behavior of the SW and SF 
networks can be explained within a mean-field picture. 
$ and T are proportional to the probability that a dis- 
placement of a random walker with diffusion constant 
C(e 2 ) is larger than 1 after D, the network diame- 
ter, time steps. The crossover size is determined from 
the condition e 2 D(N c ) ~ 1, and the spanning tree is 
very stable against the disorder as long as A -C N c . 
So, the exponential divergence of the crossover size A c 
and hence the extreme stability of the spanning tree 
are the direct consequence of the small- worldness, i.e., 
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D(N) ~ In N JyJ and D(AT) - In In In A ||, for 
the small-world and the scale-free network, respectively. 
The diameter of the SF networks with 2 < 7 < 3 scales 
as D(N) ~ In In N |2^] and within the mean- field picture 
one expects the crossover size of these networks to scale 
like N c ~ exp[exp[— ae~ 2 ]]. 

Another consequence of the disorder is that the short- 
est path from one vertex to the other becomes unique. 
Recently, it has been reported that the load distribu- 
tion in SF networks follows a power-law distribution 
with the universal load-distribution exponent 77 ~ 2.2 or 
2.0 |l(|. Phenomenologically, the SF networks with 
degenerate shortest paths seem to have 77 ~ 2.2 [To] , 
whereas those without or only a few degenerate shortest 
paths seem to have 77 ~ 2.0 0, It would be interest- 
ing to study the load distribution in the disordered SF 
networks, which will shed light on the role of the degen- 
eracy on the universality of the load distribution. This 
work is in progress. 



Acknowledgments 

We would like to thank B. Kahng for useful discussions. 
This work has been financially supported by the Deutsche 
Forschungsgemeinschaft (DFG). 



APPENDIX A: RANDOM WALK IN THE 
PRESENCE OF ABSORBING WALL 

Consider a discrete-time random walk in one dimen- 
sion in the presence of absorbing walls at x = and 
x = L. The walker is killed upon a collision with the 
absorbing walls. Initially it is at x — a and its posi- 
tion after t steps is given by x{t) = a + JZi=i x i with 
random variables Xi's with (x^ = and (x 2 ) = a 2 . In 
the time-continuum limit, the probability density P(x, i) 
satisfies a diffusion equation dP/dt = \xd 2 P jdx 2 with a 
diffusion constant /1 = er 2 /2. The absorbing walls impose 
boundary conditions P(0,t) = P(L,t) = 0, and the ini- 



tial condition implies P(x,0) — 8{x — a). In the Fourier 
series expansion P(x,t) — Y?£LiPn(t) sin^, the dif- 
fusion equation becomes dp n /dt = —(^nr 2 /L 2 )n 2 p n with 
the solution p n (t) = p n (Q)e-^ 2 ^ . The coefficients 
p n (0) are determined from the initial condition. Using 
the Fourier series expansion of the delta function, we ob- 
tain that 

2 v-^ / nira\ / nirx\ ,,/„2 t , . . 
P(x,i) = -£(sin — )(sin— je-""*. (Al) 

n=l 

The survival probability S(t) of the random walker is 
given by the spatial integral of the probability density: 

s(t) = £ sin((2n + l)7ra/£) c _ (7r 2 M/L 2 )(2ra+1) 2 t 
U 7T ^ 2n + 1 

n— 

Then the death probability during t and t + dt is given by 
S(t)-S{t+dt) ~ D(t)dt with D(t) 



dS 
dt 



Therefore the 



mean life time of the walker is given as r = J Q dt tD(t) 



f™dtS(t). Using fly), one obtains that 

8L 2 ^ sin((2n+ l)na/L) 



E 

n=0 



(2n + l) 3 



(A2) 



It becomes simpler in a special case. When a = L/2, 

r = 4gr/?(3) with (3(s) = E^oC- 1 )"^" + i)"'- Usin S 
/3(3) = tt 3 /32 



one obtains that 
T ~ 4a 2 ' 



(A3) 



When a <C L, one can use an expansion sin(n7ra/^4) 



nira/A, which yields that r 



8a L 



iL 



dd n 



Using 



that Eo 



dd n 1 



— 4C(2) = Tj- with Riemann zeta func- 



tion C(s) = 53„ =1 n s [E3l, one obtains 



ah 



(A4) 
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